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Intercluster Correlation in Seismicity 

Mega et al. proposed to use the "diffusion entropy" 
(DE) method to demonstrate that the distribution of 
time intervals between a large earthquake (the main- 
shock of a given seismic sequence) and the next one 
does not obey Poisson statistics. We have re-analyzed 
the same catalog of seismic events in California and per- 
formed synthetic tests which show that the DE is un- 
able to detect correlations between clusters, thus negat- 
ing the claimed possibility of detecting an intercluster 
correlation. We have generated synthetic earthquake cat- 
alogs with a Poisson distribution of mainshock times, 
with aftershocks triggered in direct lineage by each main- 
shock following Omori's power-law decay of seismic rate 
n(t) oc l/(t — ti) p (with p = 1, and where U is the 
mainshock time) and with a power law distribution of 
aftershock cluster sizes with exponent b/a = 1.25 (justi- 
fied from the interplay between the power-law distribu- 
tion of energies with exponent 6=1 (Gutenberg-Richter 
law) and the power-law increase of aftershock produc- 
tivity with the mainshock energy of exponent a = 0.8 
2] ) . Fig. (crosses) shows that our synthetic DE with- 
out intercluster correlation is essentially undistinguish- 
able from the real data. Thus, the conclusion of Mega et 
al. that the clusters of seismicity are correlated is not 
warranted by their analysis. However, the fact that inter- 
cluster correlation exists is known at least since Kagan 
and Jackson and can be observed straightforwardly, 
as shown in Fig. [5] 

We have used the ETAS model to generate realistic 
synthetic catalogs. The ETAS model is a simple branch- 
ing model of seismicity |^ which contains the Gutenberg- 
Richter, the Omori law, the productivity law, cascades of 
multiple triggering between earthquakes and a poissonian 
seismicity background. In the ETAS model, any earth- 
quake may trigger other earthquakes, without arbitrary 
distinction between foreshocks, aftershocks and main- 
shocks. The ETAS model reproduces many properties of 
seismicity, including realistic foreshock sequences 0. It 
is widely used to model and predict the spatio-temporal 
distribution of seismicity (e.g., Q and ref. 8 of 0]). As 
shown in Figs. 1 and 2 (circles), the ETAS model recov- 
ers (i) Mega et al.'s observation of DE S(t) = A + 5lnt 
with (5 = 0.94, (ii) the power-law pdf of the time in- 
tervals t^ 1 ™] between two successive mainshocks (which 
arises because many mainshocks are also aftershocks of 
other mainshocks) and (iii) a correlation between 's 
(which is absent in Mega et al.'s LR model). Thus, their 
LR model, which introduces correlation between clusters 
of seismicity, is insufficient to account for the correlation 
observed in the data. 

We believe that these discrepancies with Mega et al.'s 
conclusions stem for their incorrect use of the pdf of 
intracluster interevent times tj+i — tj. Omori's law is 
the pdf of the times between a mainshock and its af- 
tershocks and describes a non-stationary process such 
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□ Poisson process 5=0.5 
X GP model 5=0.92 

ETAS model 5=0.93 
DATA 5=0.94 




FIG. 1: The DE S(t) defined in Qj calculated for a Pois- 
son process (squares), the California catalog used in pj (di- 
amonds), the GP model used in and the ETAS model 
(circles) with parameters: branching ratio n — 0.93, b = 1, 
Omori exponent p = 1.3 for first-generation triggering and 
a = 0.8. 
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FIG. 2: (a) probability density function (pdf) of the time 
intervals r' m ' between two successive large (M > 5) earth- 
quakes in the Southern California catalog (circles) and in the 
ETAS model (crosses), (b) recurrence plot In r|™J vs In r 4 . 
The straight line is the linear fit to the data (circles). The 
correlation is r — 0.43 (significance > 99.9%) for the data and 
r = 0.21 (significance > 99%) for ETAS. 



that the interevent times increase as time increases since 
the mainshock: for p = 1 for instance, tj + ± — tj ~ tj. 
Therefore, the pdf of interevent times is not equivalent 
to Omori's law and it looses most of the information on 
correlations between aftershocks. This work is supported 
by NSF-EAR02-30429 and by the Southern California 
Earthquake Center (SCEC). 
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